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Abstract 

We investigate the geometry of the critical fluctuations for a gen- 
eral system undergoing a thermal second order phase transition. Adopt- 
ing a generalized effective action for the local description of the fluc- 
tuations of the order parameter at the critical point (T = Tc) we 
show that instanton-like conflgurations, corresponding to the minima 
of the effective action functional, build up clusters with fractal geom- 
etry characterizing locally the critical fluctuations. The connection 
between the corresponding (local) fractal dimension and the critical 
exponents is derived. Possible extension of the local geometry of the 
system to a global picture is also discussed. 



In a system undergoing a second order phase transition, coherent fluc- 
tuations at all scales, occur, at the critical point. The understanding of 
the geometry of these fluctuations is a long-standing problem |T]]. The self- 
similarity of the critical fluctuations suggests the formation of clusters, with 
nonvanishing value of the order parameter, which have a fractal structure 
[0, 0]. Simple geometrical arguments lead to the conclusion that the frac- 
tal dimension, describing the critical clusters, reflects the scaling properties 
of the underlying fluctuations and therefore can be related to the critical 
exponents characterizing the phase transition P]. Rigorous mathematical 
derivation and a deeper understanding of the origin of such a relation is 
however missing in the general case. Some important efforts have been per- 
formed concerning the geometry of the critical clusters in the Ising and Potts 
models ^ where the order parameter is described through a discrete field 
variable. For the more general case of continuous fields there is however, 
to our knowledge, no such understanding. It is the purpose of the present 
paper to illuminate the way the geometry of the critical clusters emerges 
in the case of a continuous effective field theory. For our considerations we 
study a self-interacting scalar field (j) at thermal equilibrium. The effective 
action of the thermal system, at the critical point T = oi the continuous 
phase transition, can be specified, in a wide class of critical phenomena, by 
an effective theory in d dimensions, in terms of a macroscopic field (j) (order 
parameter) as follows: 

r,[0] = grh-'-' j d'x ^-{V,4>f + (72A''+'|A-'^0|'+'] (1) 

The dimension of the field in (|lj) has been chosen: ~ (volume)"^ and 
the ultraviolet cut-off A of the underlying microscopic theory fixes the coarse 
graining scale Rc ~ A"^ of the effective system (throughout this work we 
use the convention = 1 (Boltzman constant) and the energy is given in 
inverse length units). The form (|I]) for constant fields leads to the standard 
equation of state at T = T^: ^ ~ (j)^ {(j) > 0) and therefore the index 6 is 
identified with the isothermal critical exponent of the system. For a high- 
temperature phase transition, the coarse graining cutoff A is bounded by the 
critical temperature itself (A ^ Tc) and the dimensionless parameters gi, g2 
in eq.(0) are expressed in terms of the ratio A = ^ (A ^ 1). In fact, using as 
a concrete example the 0{N) 3d effective theory, the action Tc[4>] in the large 
limit and for a fixed orientation in the internal 0{N) space, is written as 
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follows §: 
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It belongs to the general class (|1|) with: d = 3, S = 5, gi = X^, g2 = 2 (^7^) 
and to the particular sector gi ^ 1. 

The scalar field (j) iii the following will be no further specified. It can 
describe magnetization density or particle density or the density of any other 
extensive physical quantity characterizing the phase transition. Introducing 
the dimensionless quantities: (p = A~'^(j) ; Xi 
effective action (|l|) as: 



Axi we can rewrite the 
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In what follows we use, for simplicity, the old notation (0, Xi) instead of 



The statistics of the critical system is resolved if we are able to calculate 
the partition function: 

Z = [ P[0]e-r=[<^] (3) 



The nontrivial task is to carry out the path integration in (^. Since the self- 
similarity is a subtle symmetry of the system and it is expected to dominate 
in the formation of the fractal geometry of the clusters (the precise definition 
of the cluster will be given later on) the conventional methods to perform the 
path integration, which can be found in the literature , are not suitable for 
the problem at hand. We propose therefore to perform the summation, with 
an appropriate measure, over a class of saddle point configurations which are 
expected to dominate the critical fluctuations in eq.(0) for gi ^ 1. In order 
to illustrate our method we will treat, for simplicity, the one-dimensional case 
in some detail. Our approach can be however extended to higher dimensions 
without difficulties. 

Let us now become more quantitative. In the one-dimensional case the 
partition function of the critical system (T = T^) is given as: 
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where R is the size of the considered system. We assume here that our 
investigations refer to an open subsystem of the entire physical system located 
in the vicinity of the point a; = 0. Therefore no restrictions to the values 
of the order parameter at the boundaries of the subsystem are imposed. 
The geometrical properties of the subsystem R are expressed through the 
scaling properties of the extensive quantities characterizing the subsystem 
with varying size R around x = 0. 

The saddle-point configurations (/){x) fullfill the Euler-Lagrange equation 
corresponding to the effective action Tc[(p] and describe the classical motion 
in the concave potential U{(f)) ~ — This equation can be solved an- 
alytically in terms of two parameters E and 0(0), where ii^ is a conserved 
(during the classical motion) quantity identified with the total energy of the 
moving particle: 



Using eq.(|^) one can show that configurations with E contribute to the 
partition function Z with a suppresion factor e"^^"^'^' suggesting that the 
dominant saddle points contributing to (^ come from those solutions of the 
equations of motion for which E ^ 0. In fact eq.(^ can be integrated to 
give, for E = 0, instanton-like solutions of the form: 
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with c = 
to: 



Setting Xo = -F~7^0(O) 2 the 0-field for E = simplifies 



0(x) = A\x — Xo\ 



A 



92 



(7) 



To perform the path integration in (1) we have to sum up the contribu- 
tions of all instanton-like saddle point configurations of the form (|^), i.e. to 
integrate over the parameter Xo- In order to determine the correct integra- 
tion measure we consider the class of solutions (H) with Xo 3> R. In this 

2_ 

case (j){x) =constant ~ Xo and in a region of radius R around x = 

the path integration becomes an ordinary integral over Xo with measure: 
4+1 

Vcp = djj,{xo) ~ Xo dxo- To determine the range of integration over Xo 
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we have first to clarify tlie meaning of a cluster in our picture: Let us as- 
sume that M is an extensive variable (i.e. magnetization) characterizing the 
field configurations of the critical system and possesing a minimal, in general 
different from zero, value (M > //) related to the microscopic details of 
the system. Within the picture of the local observer positioned at x = a 
cluster of size R is the set S of points with a maximum distance R from the 
origin. The most appropriate observable to study the geometric properties 
of the critical clusters, is the thermal average < M{R) >=< (p{x)dx > 
and in particular its behaviour as a function of R. The minimum value /i 
of the magnetization M introduces a threshold (pmin to the configurations 
contributing to this average, leading us to an upper limit for the integral Xo 

S-l 

as a function of R: Xo < (^) ^ • Oii the other hand the singularity at 
X = Xo must lie outside the region of the considered cluster (0, R) restricting 
the integration in (§) over configurations with Xo > R- In terms of the 
partition function (|^) this average can be determined as: 

fR AAR/p,)'-^ _s+i ,_3 

< / (j){x)dx > = / dxo ■ Xo \xt' - {xo - R)~] (8) 

R 

^[(Xo-R)'— - Xo 

with Gi = 2gig2A^+\ O ne can show analytically that in the large Gi limit 
{Gi ^ 1) there are three characteristic regions determining the behaviour of 

S+l S+1 - 1^-1 

the integral in eq.(8). Putting R^ = A~~s~ fi~s~Gi and Ru = Gf"*"^ we find 

that in the central region of scales <^ i? -C Ru'- < (p{x)dx >~ 
leading to a fractal structure of the cluster around the point x = with a 
fractal mass dimension 0: 

^^ = 5TT 

This behaviour crosses over for i? 3> -Ru to a different power-law: 
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rR 

< / (j){x)dx >~ 



suggesting the presence of a fractal with mass dimension dp = jzj at large 
scales. 
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For R <^ Rd a. violation of the scaling symmetry of the critical cluster 
is revealed leading to an approximately constant value of the integral (8). 
The parameter R^ defines a minimal scale of the critical system effectively 
related with the minimal value of the order parameter. In Fig. la we show 
the numerical results for the calculation of (8) using the values Gi = 5 • 10^ 
and 5 = 5. The three different regions and the corresponding cross-over scales 
describing the geometry of the critical cluster are clearly distinguished. In 
the same plot we show also the corresponding linear fits to illustrate more 
transparently the above considerations. The fractality in the central region 
characterizes the critical system in the sense that it corresponds to the scaling 
behaviour in the vicinity of the local observer when ;U — > 0. This is at best 
shown in Fig. lb where we calculated (8) using the same value of Gi as in 
Fig. la and let the upper limit in the Xq- integration going to infinity. The 
crossover scale R^ gives presumably a measure of the correlation length of 
the finite system at T = T^. 

One can now easily generalize the above investigations in order to describe 
systems of higher dimensions [d > 1). We proceed in a similar way as for the 
one-dimensional case using the saddle point approximation for the partition 
function Z^. 

Zd = j P[0]e-^^i'v"'"'[^(^^'^)'+^^l*l'^'l (9) 

in order to calculate the thermal average < /q^ dl^x >Za- The summation 
over the saddle points in eq.(^ becomes an ordinary integration over the 
position of the singularity in the instanton-like solutions <pd in ^-dimensions. 
In an analogous way as in the Id case, solving the Euler-Lagrange equations 
for the critical action, we get instanton-like solutions for (pd- In the case d = 2 
a class of analytic solutions possessing a point singularity can be determined: 

-1 

d=2 ■ ^,(f^ = A,\r-r,\^^ ; = (|(5 - 1)^)'" (10) 

Performing the calculation of the mean value : < M{R) >=< J d'^f(j){'f) >, 
characterizing a two-dimensional critical cluster, in an analogous way as 
for the Id case, we get a similar behaviour concerning its fractal geomet- 
ric properties. There are characteristic scales in the radial component Rd = 

A2^fi^Gf and i?„ = with G2 = Tr9i[j0^ + such that: 

d = 2 ; < M{R) > ~ R^^ ; Rd <^ R <t: Ru 
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< M{R) > ~ ; Ru-^R (11) 

A cross over for large R is found also in this case. For dimensions d > 3 no 
analytic solution to the Euler-Lagrange equations is available in the general 
case of a non vanishing anomalous dimension r] (5 = [^). However 

these equations can be integrated numerically leading again to an instanton- 
like behaviour. In particular one can find exact analytic spherical solutions 
of this kind for > 3 in the special case when rj = {S = 

For < 77 ^ 1 an approximate solution can be obtained given as follows: 



d-2 

2-d f {d — 2)ro\ 2 




(12) 

The solution (|12D goes to the exact one for 77 = 0. For r far from the 
singularity region the approximate form ( |T^ ) coincides practically with the 
exact (numerically obtained) solution. This can be at best seen in Fig. 2a 
where we plot together the numerical and the approximate solution to the 
Euler-Lagrange equations for d = 3 and r] = 0.34. In fact, in a wide range 
of universality classes including the 0(4) theory in which rj ^ 0.034 [0], 



the anomalous dimension for d = 3 is much smaller and therefore one 



can safely use the solution (|I2|) for most calculations. Based on (|12D we 
determined < M{R) > for spherically symmetric clusters in d > 3 dimensions 
and explored the geometric properties of such a cluster. Once again we 
got the typical central fractality region crossing over to a fractal with a 
smaller dimension for distances comparable to the correlation length. It 
must be noted that for d > 3 the cross over disapears as r] 0. Using 

Ra = a^^H^ (Gd)^ and Ru = (Gd)^^^^ with a = ' ,Gd = 

dT(d/2) 9i92 and q = — we obtam: 

d>3 ; < Md{R) > ~ R^^^ ; i?d < i? < i?„ 

< Md{R) > ~ ; K < ^ (13) 



The characteristic behaviour of < M{R) > for (i = 3 is presented in Fig.2b. 
Here we used G3 = 10^ and 77 = 0.34 (as in Fig.2a). The breaking of the 
fractality (for R -C Rd) is clearly reproduced while the cross over is sup- 
pressed due to the small value of i]. The power-laws < Md{R) >~ -R*^^ or 
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~ i?*^^ with d = 1,2,.. determine fractals at different scales with dimensions 
dp, dp- Putting together our results in 1,2 and 3 dimensions and taking into 
account that our considerations for d > 3 are restricted to the case when the 
anomalous dimension rj is small (?7 ^ 1) we can cast our results for the frac- 
tal properties of the critical cluster into universal expressions determining dp 
and dp in terms of d and 6: 

"--—1 • = 1 • 2(2-,) 

The expression for dp is in accordance with the results obtained in [0, Q 
for the Ising and Pott's critical clusters. Thus we have found that eqs. (|I^) 
describe the fractal geometry of the critical clusters in a wide range of scales 
and for a general class of effective theories. Considering the asymptotic re- 
gion where the size of the cluster reaches the value of the correlation length 
of the finite system we find a cross over to a more dilute phase with a smaller 
fractal dimension dp. Furthermore we have revealed a mechanism responsible 
for the formation of the geometry of the critical clusters. Our considerations 
are restricted to the point of view of a local description. The generalization 
of our approach in order to build up the entire critical system, requires the 
extension of the formalism to configurations incorporating many, suitably lo- 
cated, instanton-like structures (of the size of the correlation length) covering 
the whole available space. Such an investigation however goes beyond the 
scope of the present letter. 
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Figure captions 

Figure 1: (a) The mean magnetization < M{R) > as a function of R (in units 
A~^) for the d = 1 case, Gi = 5 ■ 10^ and A = 1. The fitted lines indicate the 
two regions of fractality as described in the text. A nonvanishing minimum 
magnetization fi takes care for the violation of the scaling in small distances 
R. 

(b) The mean magnetization < M{R) > for d = 1, Gi = 5 ■ 10^ and A = 1 
for /i ^ 0. 

Figure 2: (a) A saddle point solution to the action (|1]) for c? = 3 and t] = 0.34. 
The analytical approximation (dashed line) and the result of the nummerical 
integration (solid hne) are displayed separately. 

(b) The mean magnetization < M{R) > for d = 3 and = 10^. The fitted 
line indicates the fractality in the central region. 



9 



(a) 




4n 



(a) 



3- 



numerical solution for ri=0.34 
analytic approximation 



2- 



1 - 



0,0 



0,5 



1,0 



1,5 



2,0 



(b) 

10' -j 



106 -J + <M(R)> for ri=0.34 

linear fit with slope s=2.23 ^ 

105 i 



100 = 



104 = 



103 = 



10^ ^ 



101 ^ 



Figure 2 



10 100 1000 



R 



